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It is shown that the power spectrum defined in the Synchronous Gauge 



^S| ■ can not be directly used to calculate the predictions of cosmological models 

^N , on the large-scale structure of universe, which should be calculated directly 

r — . by a suitable gauge-invariant power spectrum or the power spectrum defined 

o. 

in the Newtonian Gauge. 
O 

:^ 
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^ : PACS: 98.80.-k; 98.80.Hw; 98.80. Cq 



It is indispensable to adopt the linearized Einstein equation to solve the problem of the 
evolution of density perturbation in the expanding universe, and this involves the freedom 
of gauge. Currently, the density perturbation is widely calculated in the Synchronous Gauge 
{SG) frame and then directly compared with the observations of the large-scale structure 
of universe, e.g., the abundance and distribution of the galaxies or the clusters (e.g., Refs. 
|jl|-§|]). Because the power spectrums of density perturbation differ dramatically for small 
wavenumbers in different gauges, this raises an important problem: Does the density per- 
turbation in the SG frame characterize directly the large-scale structure of universe which 
is relative to the unperturbed expanding background of universe? I do not think so. 

Because the density perturbations defined in different gauges are generally different from 
each other, we need to construct some suitable gauge- invariant variables that can be directly 
confronted with the large-scale structure of universe. The gauge-invariant approach to the 
gravitational perturbation, which was pioneered by Bardeen [H, has been applied to the 
cosmological perturbation (e.g., Refs. |0J^). However, there exist some errors in the com- 
parison between the theoretical calculations and the observations of the large-scale structure 
of universe. 

The most general form of the line element for a spatially flat background and scalar 
metric perturbations can be written as [§ 

ds"^ = a^Kl + 2^)dT^ - 2B\idx'dT - [5,j - 2{<^5ij - EQ]dx'dx^ }, (1) 

where a and r are the conformal cosmic expansion scale factor and the conformal cosmic 
time; "[" denotes the background three-dimensional covariant derivative. The corresponding 
perturbed energy-momentum tensor T^ has the form 

T\ = p{l+e), 
T,' = {p + p)U\,, 

TV = -p(l+£W.-(p + p)S|,, , (2) 



here p and p are unperturbed density and pressure; U and E determine velocity perturbation 
and anisotropic shear perturbation, respectively; e and e denote the fluctuation for density 
and pressure respectively. 

In many gauge-invariant approaches, the gauge-invariant density perturbation was con- 
structed by @,|^ 

em = e + ^{l + w)-{U-B) , (3) 

a 

where a dot denotes the derivative with respect to the conformal time t; w = p/p is the ratio 
of the pressure to the density of background, e^ coincides with the density perturbation 
^(CTG) in the Comoving Time-orthogonal Gauge {CTG, in which U = B = 0), which denotes 
the density perturbation relative to the spacelike hypersurface which represents the matter 
local rest frame everywhere ^. This quantity also coincides with the density perturbation 
€(^SG) ill the Synchronous Gauge (in which \I^ = i? = 0) for the pressureless matter system. 
In other words, e^ denotes the density perturbation relative to the observers everywhere 
comoving with the matter {not with the unperturbed expanding background of universe 
because the matter have local velocity other than the Hubble flow). So, e^ has physical 
significance only for the small-scale perturbation everywhere, and this can be also derived 
from following facts: The relation between e^ and the general gravitational potential (p 
(which was correctly constructed in literatures, see Eq.(^ below) obeys Poisson equation 

m-- 

y2 ^ _/.20 ^ AnGpa^em ■ (4) 

where k is the (comoving) wavenumber of Fourier mode. Poisson equation is valid only 
in the case that the investigated scale is small compared to the Hubble radius 1/H with 
H = d/a^ being Hubble constant Q |]T|. Hence, e^ can not be regarded as the physical density 
perturbation directly for the large-scale (comparable to or larger than the Hubble radius) 



^In the standard big-bang model, the Hubble radius is the same order as the horizon, but in the 



perturbation mode in the expanding universe because the gravitational potential has been 
constructed correctly. Another fact is that the linearized Einstein equation shows that e^ 
keeps on growing even for the perturbation modes outside both the Hubble radius and the 
horizon ||^,^ , which also means that e^ can not directly characterize the physical large-scale 
perturbation modes. Therefore, the physical density perturbation which comprises both the 
small-scale and the large-scale modes can not be described directly by Cm, though it is a 
gauge-invariant quantity. 

In fact, the large-scale density perturbation has physical significance only in the case that 
the corresponding density perturbation everywhere is relative to the unperturbed expanding 
background of universe. So we must construct a proper gauge-invariant variable to describe 
the physical density perturbation directly^ and then confront it with the large-scale structure 
of universe. It is my viewpoint that the gauge-invariant quantity constructed by 



eg = t 



^(B-E) = e^3(l + w)-(B-E) . (5) 

p a 

is just what we sought for. eg coincides with the density perturbation e(iVG') in the Newtonian 
Gauge [NG, in which B = E = 0). Correspondingly, two gauge-invariant scalar potentials 
and ip, both of which become the same as the gravitational potential in the Newtonian 
limit, are constructed from metric perturbations |^: 

^ = <^--{B-E) , 
a 

^^^i; + ]L^[^B-E)a]. (6) 

a cLt 



inflationary scenario, the horizon can be much large than the Hubble radius. The largest scale in 
causality is set by the horizon, while the physical perturbation modes outside the Hubble radius will 
hardly grow as long as the equation of state for the background of universe is unchanged. Notice 
that the perturbation amplitude for all modes will change if the equation of state of the background 
changes, e.g., (p changes during the transition from radiation-domonated era to matter-dominated 
eraH. 



The gauge-invariant pressure perturbation 5p/p and velocity perturbation Vi are constructed 

as 

5p/p^e + ^-{B-E), (7) 

P 

v, = {U-E)\,. (8) 

The time-time part of the hnearized Einstein equation gives [^,^ 

y2 _ 3^(^^ + 0) = _A;20 _ 3^(^^ + 0) = 47rGpa2e„ . (9) 

a a a a 

From this equation we can see that the relation between the potential (p and the physical 
density perturbation eg does not obey Poisson equation any more for the large-scale pertur- 
bation modes. At the same time, it can be deduced from the linearized Einstein equation 
that eg hardly grows when the perturbation mode is outside the Hubble radius {k < aH). 

Hence, though any gauge can be employed to work for the density perturbation. We 
should adopt the gauge-invariant variable which characterizes the density perturbation rel- 
ative to the unperturbed expanding background of universe to compare with the large-scale 
structure of universe directly (here I do not consider the biasing issue). 

For the mult i- component system, e^, 5p/p and f , denote the gauge- invariant variables of 
total density perturbation, total pressure perturbation and total velocity perturbation re- 
spectively, and the corresponding gauge-invariant variables for a-component are constructed 
in similar ways [0: 

e^g = eo. + ^{B-E) , (10) 

Pa 

{Sp/pU^e^ + ^{B-E), (11) 

Pa 
V^, = {Ua - E)\i . (12) 

From above equations we can see that all these gauge-invariant perturbation variables 
conincide with the corresponding ones in the Newtonian Gauge respectively. We can also 



adopt the density perturbation defined in this gauge to compare with the large-scale structure 
of universe directly. 

Any a real-space fiuctuation 5(x, r) can be decomposed into Fourier modes 

5(x,r) = -i-^/4(r)e^'^-^rf^fc, (13) 

[271)2 J 

with the reality condition 61- = (5_k . The power spectrum is defined as the mean square of 
the corresponding Fourier-mode amplitude: 

Ps{k,r)=<Mr)\'>^\6{k,T)\\ (14) 

All the statistical properties of the Gaussian random field can be determined by the power 
spectrum completely, and for simplicity we only consider this case here. At the same time, 
we are only interested in the fiuctuation relative to the unpertured background everywhere 
which can be confronted with the large-scale structure of universe directly. The root mean 
square (rms) fiuctuation and two-point correlation function of density on the scale R can be 
calculated by the physical power spectrum P^ (k,T) = |eg(A;,r)p directly 



1 /-oo fjU 

^..(^,r)=[— /^ P.^{k,r)W\kR)e- 

sin kR , q dk 



^ /^ N I f^ ^ ,, .sin kR,ndk ,,^. 

where, W{x) = -ji{x) is the window function which truncates the contribution of the 
power spectrum of large wavenumbers to the rms fiuctuation (here ji is the first-order 
spherical Bessel function). When calculating <Jeg{R, t) or C,eg{R, t), we need to integrate the 
power spectrum over all wavenumbers, not merely the modes within the Hubble radius (see 
Eq.(|T5])). Because the spectrum P^^{k,T) = \em{k,T)\^ differs dramatically from P^^^k^r) 
for small wavenumbers (e.g., k < aH), generally we can not replace P^g{k,r) with P^^{k,T) 



in Eq. (|1^) directly to obtain the rms fluctuation and correlation function of the density. 

The temperature fluctuation of the cosmic microwave background (CMB) can be char- 
acterized by the angular power spectrum C/. Based on the line-of-sight-integral method 



developed by Seljak and Zaldarriaga [10], we have 
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6(A;, T)ji{kTo — kT)dT 





(16) 



with 



e(A;, r) = g (^e, + 0+| + ^|j+ e-«(0 + ^) 

where j/ is the spherical Bessel function; e^ is the photon density perturbation; H is a gauge- 
invariant quantity denoting the effect of polarization and anisotropic stress of photons; Vb 
denotes the velocity perturbation of baryons; tq denotes present time; g is the visibility 
function and k, is the total optical depth at time r. 

Now, let us take the standard cold dark matter (SCDM) model as an example to illustrate 
the prediction differences between the conventional treatment and my viewpoint. The SCDM 
model (the cold dark matter density parameter Q^ = 0.95, the baryon density parameter 
Qf, = 0.05, Hq = lOOh km s~^ Mpc~^ with h = 0.5) assumes that the primordial density 
spectrum in the SG frame Pe{SG){k,Ti) is the Harrizon-Zel'dovich (HZ) spectrum: 



<SG) 



{k,r,)=\e^SG)ik,Ti)\^=Ak, (18) 



here A is a constant and Tj is the primordial time. In the conventional treatment, the 
density perturbation ^(sg) iii fhe SG frame was directly regarded as the physical density 
perturbation, so the rms fluctuation of the density on the scale of rg = 8h~^ Mpc was 
directly corresponded to (Te(5G)(^8, ^j), which is related to Pe{SG){k,To) =\ e(sG){k,To) p by 



(^e{SG)ir8,ro) 



1 /"^^ Hk 

^/^ P.^SG){k,ro)W^krs)e- 



(19) 



2 

In the foregoing discussions we have shown that this picture is wrong because e(^sG){k, t) 
(or its spectrum P^i^sG){k ■, t) =| e(^sG){k,T) p) can not characterize the physical density 
fluctuation directly. The physical density spectrum corresponding to the assumption of the 
SCDM model on the primordial density spectrum should be 

Pe^k^r) = P.^{k,n)Tl{k,T) ^ [M^] AkT^^ik,T) , (20) 



where, the subscript 'V' denotes the primordial time; P^^{k,Ti) ~ {3ajHf/k'^)^Ak is the 
physical primordial spectrum corresponding to the HZ spectrum in the SG frame (notice 
that all the interesting modes are far outside the Hubble radius at time Tj, i.e., k << aiHi), 
and T^g{k,T) = ^^^^ '^l is the transfer function of physical density spectrum due to the 
evolution of density perturbation from the primordial time Tj to the time r. 

Fig.l shows the physical density spectrum P^g{k,T) of the SCDM model at the redshift 
z = and z = 10, with the corresponding P^(^sG){k,T). It can be seen that Pe{SG){k,T) can 
approximate the physical spectrum only for the large wavenumbers (small-scale modes), and 
differs significantly from the latter for the wavenumbers outside the Hubble radius. By the 
way, Pe{SG){k,T) differs hardly from the gauge-invariant power spectrum P^^{k,T) for the 
SCDM model in which the universe was dominated by the pressureless matter at very high 
redshift. In the conventional treatment, only the power spectrum of the wavenumbers inside 
the Hubble radius is taken seriously, and the part outside the Hubble radius is thought to 
have no physical significance. This viewpoint is not correct: When a real-space fluctuation 
is decomposed into Fourier modes, all these modes will have physical significance whether 
or not they are inside the Hubble radius (see Eq.(|T3D). Only the whole power spectrum can 
determine a real-space Gaussian fluctuation completely. In fact, the power spectrum of the 
wavenumbers outside the Hubble radius has also been included in the current calculation of 
the density fluctuations (e.g., see Eq.(|l9])) and the angular power spectrum Ci for CMB (see 
Eg . (p!6D ) . So, the conventional treatment is inconsistent upon this point. More important, 
the temperature fluctuation detected by the COBE satellite (e.g., C2) is dominated by the 
power spectrum of the wavenumbers outside the Hubble radius at the recombination era! 
On the other hand, inflationary scenario has provided a well-known mechanism of the non- 
vanishing power spectrum for the wavenumbers outside the Hubble radius. 

Fig.l also reveals another important difference between the physical spectrum and 
Pe{SG){k,T): the shape of the physical spectrum changes even for the pressureless mat- 
ter system, because the physical density perturbation modes far outside the Hubble radius 
can not grow due to the causality; while the shape of Pe(SG){k,T) for the SCDM model 

8 



hardly changes, because the growth rate of e(^sG){k,T) is independent of wavenumber k for 
the presureless system. 



Now we re-check some predictions of the SCDM modeL From Eqs.(|l5|) and (|20| ) we 
find that the rms fluctuation a^g {R, t) and the correlation ^^^ (-R, t) of density are both 
divergent on any scale R at any time r, so do those of the potential 0, if the physical density 
spectrum has the HZ shape exactly in the limit oi k ^ in the SG frame. In fact, these 
unreasonable predictions caused by the divergent spectrum have also been reflected in the 
calculations of the temperature fluctuation of CMB, though it is very implicit. From the 
combination of Eqs. ([T6|) , (0) and (^0]) we can obtain that Cq is divergent too. This means 
that the temperature fluctuation of CMB is infinite in any direction. Because people are 
only interested in the angular distribution of the temperature anisotropy of CMB (e.g., Ci 
for / > 1), this unreasonable prediction of the SCDM model has been ignored so far. 

All these unreasonable predictions are due to the fact that the power spectrum was not 
cut off effectively. Any physical power spectrum should become negligible as wavenumber k 
tends to zero due to the causality. So, the primordial spectrum should be cut off effectively 
at a minimum wavenumber kmin, e-g-, CL/kmin is at most the same order as the horizon (notice 
that the horizon may be much larger than the Hubble radius). After doing so, we can see 
that all calculation results are convergent and dependent on kmin- 

The key point lies in how to determine the effective minimum wavenumber fc^m for the 
power spectrum. If kmin is inside the Hubble radius at any time, there will not be so many 
problems to discuss here. However, kmin can be far outside the Hubble radius (for example, 
the quantum fluctuation in the early universe can be stretched far outside the Hubble radius 
by inflation process), otherwise we can not explain the temperature fluctuation detected 
by the COBE satellite, which showed that these fluctuation modes are outside the Hubble 
radius at the recombination era. 

In conclusion, the power spectrums differ dramatically for the modes outside the Hubble 
radius in different gauges, however, the calculation of the rms density fluctuation needs 
to integrate the power spectrum for all perturbation modes (see Eq.(^)), not merely the 

9 



modes within the Hubble radius, so do the calculations of the two-point correlation function, 
the abundance of the galaxies or clusters by Press-Schechter formula, and the angular power 
spectrum Ci for CMB, etc. In the above demonstrations we have shown that these quantities 
should be calculated directly by the gauge-invariant spectrum P^g{k,T), or the spectrum 
Pe{NG){k,T) =1 e(^NG){k,T) p in the Newtonian Gauge which coincides with P^^(k,T), but 
not directly by the spectrum P^(^sG){k, r) in the Synchronous Gauge. On the other hand, the 
density fluctuations and the CMB anisotropies are dependent on the cut-off wavenumber 
kmin of the physical power spectrum, kmin is not set by the Hubble radius simply, but by 
the physical process in early universe, such as inflation. In the inflationary picture there 
generally exists a non-vanishing power spectrum for the modes outside the Hubble radius, 
i.e., kmin can be far outside the Hubble radius! These conclusions will change dramatically 
the scenarios of the current cosmological models on the formation of the large-scale structure 
of universe. 
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Captions for Figure 

Fig.l The physical density spectrum for the SCDM model (denoted as the solid lines), 
with the corresponding spectrum in the Synchronous Gauge (denoted as the dashed lines). 
Since the physical rms fluctuation cr^^ir^, tq) (see Eq.(|l^)) is divergent for the SCDM model, 
for simplicity the normalization for the density spectrum is chosen so that a^[SG)ij8, tq) = 1 
(see Eq.(0)) instead; The upper and the lower lines correspond to the cases of ^ = and 
z = 10, respectively. 
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